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Abstract 

High energy ion beams (> MeV) generated by intense laser pulses promise to be 
viable alternatives to conventional ion beam sources due to their unique prop- 
erties such as high charge [U |2] , low emittance [3l [4] , compactness and ease of 
beam delivery [5 . Typically the acceleration is due to the rapid expansion of 
a laser heated solid foil, but this usually leads to ion beams with large energy 
spread. Until now, control of the energy spread has only been achieved at the 
expense of reduced charge and increased complexity [H [71 18] . Radiation pres- 
sure acceleration (RPA) provides an alternative route to producing laser-driven 
monoenergetic ion beams [9l|T0j. In this paper, we show the interaction of an 
intense infrared laser with a gaseous hydrogen target can produce proton spectra 
of small energy spread (a ~ 4%), and low background. The scaling of proton 
energy with the ratio of intensity over density (I/n) indicates that the acceler- 
ation is due to the shock generated by radiation-pressure driven hole-boring of 
the critical surface [HI [12] . These are the first high contrast mononenergetic 
beams that have been theorised from RPA [9l HOl [131 |14l [15] , and makes them 
highly desirable for numerous ion beam applications. 
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Most investigations of ion acceleration by intense lasers have relied on sheath 
acceleration. When a high-intensity laser impacts on a solid target, an over- 
dense (opaque) plasma is formed, i.e. Ue > ricrit = eomc^o/e^. Laser energy 
is converted into hot electrons with temperature ksTe^ which expand to form 
an electrostatic sheath around the target with fields of order ~ MV//im. This 
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sheath accelerates surface ions to energies ~ ksTe (~ multi-MeV) per nucleon 
[161 [17]. However, this usuahy produces ion beams of broad energy spread, with 
a large number of ion species, many originating from impurities [HI [2]. Also, 
the production of proton beams has until now mostly depended on the presence 
of the same target impurities (hydrocarbons and moisture). Due to the bulk 
heating, ions are accelerated from both front and back surfaces, providing in- 
efficient energy coupling between the laser and the ions, and ion energy scales 
only weakly with irradiance (~ (/A^)-^/^) [191 [20]. Modulated ion energy spec- 
tra have been produced, either by using complex target preparation [7l|6], or 
energy selection , but these methods still feature relatively poor peak-to-noise 
spectral contrast, reduced yields and the presence of impurity species. 

Radiation pressure acceleration (RPA) has been proposed as an alternative 
method of ion acceleration at ultrahigh intensities [9l O [El ttH [Ej • For an 
opaque plasma (n > ncr), the laser radiation pressure, Pr = 2//c, initially 
pushes plasma electrons into the target setting up an electrostatic shock whose 
space charge field pulls along ions at the hole-boring velocity [11] IT^, v^b ^ 
(2///?c)^/^, where / is the intensity and p is the mass density. Stationary ions 
in advance of the shock front can be accelerated by the same space charge field, 
effectively bouncing off the shock front, producing a population of energetic ions 
with velocity 2vhb [lH |22l |23J . If the target becomes sufficiently thin that all of 
the electrons can be pushed out of the target {d < c/ujp)^ the ions can be pulled 
along in unison, in what is called the "light-sail" phase of RPA. The "light- 
sail" phase is necessary to reach higher energy [15] |14^ . First experiments have 
now reported limited enhancement of ion energy contrast with optical lasers 
and ultrathin solid targets [24]. However, due to the extremely high intensities 
required, energy gains were modest, energy spreads large, and impurity- free 
proton beams could not be produced. 

Gas targets have been shown to be an alternative to solid foils. They operate 
at high repetition rate, are easily adjusted for target density and material and 
are less prone to contamination [25]. However, they are difficult to operate at 
near or above critical density, which is necessary for efficient ion acceleration 
(e.g. for 1 /im laser, Ucrit — 1 x 10^^ cm~^) [26^,^27|. This can be remedied by use 
of longer wavelength (infrared) lasers. For example, for a A ~ 10 /im CO2 laser, 
^crit {o^ A~^) ^ 10^^ cm~^. This density is easily obtained by ionisation of gas 
targets. Importantly, due to the lower specific mass (< 10^ x solid densities), 
gas targets become of interest for RPA at proportionally reduced laser intensity. 

An experiment designed to investigate radiation driven acceleration of gas 
targets by an infrared (A = 10 /im) laser at intensities up to 10^^ Wcm~^ was 
performed (see methods). The interaction was diagnosed by transverse probing 
and measurements of the forward accelerated ion beam. A gas nozzle of diameter 
L = 1 mm was employed, which though only L ^ lOOA, is not thin enough for 
"light-sail" RPA. The normalised vector potential was ao ~ 0.5 (see methods). 
Since the ponderomotive potential Up = a^rrieC^ ^ for this ao one would expect to 
produce hot electrons, which would promote sheath acceleration. However, this 
was mitigated to some extent by use of circular polarisation in the experiment 
[2i[l0]. 
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Figure 1: Ion spectra: Selection of raw and processed proton spectra for varying peak density 
(n) and vacuum intensity (/) showing scaling of peak proton energy Emax oq I/n. Parameter 
I/n shown to the right of the respective raw images. Shots taken with; a) / = 7.5, n = G.lricn 
h) I = 6.5, n = G.lricr, c) I = 6.9, n = 7.6ncr, d) / = 6.7, n = S.Oricr (intensities in 
units lO"*^^ Wcm~^). e) Background subtracted (solid lines) and also corrected (dashed lines) 
spectra. Heights of corrected spectra adjusted to match those of raw line-outs. Line-out 
corresponding to b) reduced 4x to fit on same scale. 



Figure 1 shows a selection of ion beams recorded with the magnetic spec- 
trometer for a range of densities and intensities. Ion beams with a peak energy 
up to Emax = 1-1 MeV were observed. The most striking observation is the 
narrow spectral width of the ion feature, especially at higher I/n. The peak-to- 
noise contrast of the spectrum is greater than 100. In figure la, the ion beam 
image is a circle comparable in size to the projected aperture size at the image 
plane, indicating that this feature is dominated by the instrument function of 
the spectrometer. Spectra (figure le), taken after deconvolving the instrument 
response (see methods), show that the narrowest observed feature (figure la) 
corresponds to an rms energy spread of only a = 4.2 %. 

For the shot peaking at 0.6 MeV (figure lb), the maximum of the spectrum 
is ~ 3 X 10^^ protons/Me V/sr, which is up to 1000 times greater than for 
previous modulated ion beams from laser plasma interactions [6l [71 [24] . From 
the transverse spread of the 4 x 10^ protons passing through the spectrometer 
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Figure 2: Proton energy scaling: as a function of expected energy due to hole-boring, 
E = 41 /nc. For Sricr > n > 4:ncr, peak proton signal increases with increasing I/n. For 
n < Aricr beams were produced with larger variation in energy (below the observed scaling), 
and are not included. No beams observed for n > Sricr- Vertical bars represent rms energy 
spread of each shot. 



pinhole, the geometrical emittance of the beam was determined to be e = 0.16 
mm-mrad, corresponding to a normalised emittance of = P^e = 8 nm-rad. 

Also apparent is that the peak energy increases with increasing I/n over 
this density range. This is shown more clearly in figure 2, which plots measured 
peak energy versus expected energy due to shock acceleration, E = ^m{2vhb)'^ = 
AI /nc. The data shows a trend of increasing energy for increasing I/n. However, 
the linear best fit implies slu I/n which is 9 times higher than expected by 
taking the vacuum focused intensity, and the peak measured densities. This 
discrepancy will be addressed further below. 

Simultaneous transverse optical probing was also performed. Although opaque 
to the infrared driver, at these densities the plasma is transparent at optical 
wavelengths. By contrast, probing the surface of a solid-density interaction 
would be impossible due to severe refraction. Figure 3a shows a shadowgram of 
an interaction for peak density n = Ancr-, at t ~ 30 ps after arrival of the main 
CO2 pulse, shortly after the end of the interaction with the second pulse. A 
large cavity has been created by the laser inside the gas target. Interferometry 
(figure 3b) shows that the plasma within the cavity is at much lower density 
(< lO^^cm"^), whilst the walls of the cavity are just above critical density. 
Hence this density-discontinuity (shock) is associated with hole-boring. The 
laser has travelled up to the critical surface, at which point it is mostly re- 
fiected. The resultant radiation pressure causes the critical surface to be driven 
inwards. Probing at later times shows the shock front moving further into the 
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target, though at reduced velocity due to energy dissipation. Features are seen 
inside the overdense region of the interaction due to the hot electrons created 
at the shock front. At later times, it is possible to see a slower evolution of the 
back surface of the gas target, which was presumably not detected by the ion 
diagnostics since the associated energy is much lower. 

To model the interaction, a series of 2D particle-in-cell simulations were 
performed. A hole-boring front can be seen (figure 4a) forming above critical 
density, but well below the peak density. The velocity of the front increases up 
to the peak of the laser pulse where it is comparable to the expected hole-boring 
velocity, v^b- As the peak density is varied in the simulations, the velocity of the 
hole-boring front is found to decrease. This is because the shock, though starting 
near critical density, reaches higher density as it moves forward. Secondary 
pulses enhance the shock structure and snowplough a greater number of protons 
to Vhb, allowing protons below the hole-boring velocity to "catch-up" with the 
shock front. This can serve to improve the fidelity of the shock structure. 

The simulations also show self-focusing of the laser pulses (by a factor > 2 in 
intensity) [29 . This in conjunction with the lower than peak interaction density 
explains the previously noted discrepancy between measured and predicted ion 
energies in the experiment. Self-focusing would be accentuated in three dimen- 
sions. To account for this a further simulation was performed with ao = 0.9 and 
shorter density scale-length. Figure 4b shows the spectrum observed in a simu- 
lated spectrometer, which models the experimental one. A clear peak is obseved 
at high energy with similar energy to that observed experimentally. Assuming 
a real beam divergence comparable to that of the simulated high energy proton 
bunch (~ 4°) would imply a total number of accelerated protons in the bunch of 
up to 5 X 10^ {c^ 0.8 nC). Such a high-current low energy spread beam of protons 
would obviously be of interest for applications as diverse as hadrontherapy [5] 
and as an injector for high-energy particle accelerators. 
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Figure 3: Optical probing: (a) Shadowgraphy of the interaction for n = 4ncr, 30ps after 
the first main pulse interacts with the gas. Laser entered from left, with the silhouette of the 
gas nozzle shown above. White dotted lines approximate initial boundaries of the gas jet. (b) 
Density profiles obtained from simultaneous interferogram 
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Figure 4: PIC simulation: (a) Proton density map and (b) proton spectra as seen by 
simulated spectrometer taken 12ps after the start of the interaction. Simulation performed 
with code OSIRIS, using 20000 x 3000 cells to simulate a 600 x 500 /im (~ 60 x 50 A) box 
with laser incident on fully ionised hydrogen plasma with density profile increasing from to 
7.5nc in 100 /im. Parameters chosen to replicate the experiment, except with a higher ao to 
account for self focusing (ao = 0.9, tx, = 8ps = 229 laser periods, circular polarisation). Both 
formation of a distinct hole-boring driven shock and an enhanced population of protons at 
V = 2vhb are observed. 
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Methods 



0.1. Laser parameters 

The experiment used the A = 10 /im, 0.5 TW peak power, circularly polarised 
CO2 laser at the Accelerator Test Facility, Brookhaven National Laboratory 
The short pulse is achieved in a Kerr cell filled with optically active CS2 liquid 
via fast polarisation switching by a 5 ps long, frequency- doubled (A = 532 nm) 
NdiYAG beam. Spectral modulation of a picosecond A = 10 /im pulse in CO2 
gas amplifiers results in splitting of the output pulse into a train of 5 ps pulses 
with 25 ps period with most of the energy contained in two leading pulses. Shots 
were taken with integrated energy in the range 2.5 - 2.9 J and focused by an 
f/3 off-axis parabolic mirror to a spot size of wq ~ 70 /i. This gives vacuum 
target intensities of 6.5 — 7.5 x 10^^ Wcm~^, or normalised vector-potential ag ^ 

eE 

0.51 — 0.55, where for circular polarisation gq = :^ 0.60 (/A^)^/^ (for / in 

mujc 

10^^ Wcm"^ and A in /im). 

0.2. Experimental set-up 

The laser was focused onto the front surface of a hydrogen gas jet from a 
L = 1 mm circular nozzle. A r = 10 ps, 2lo (A = 532 nm) NdiYAG laser beam, 
absolutely synchronised with the CO2 beam, was used for probing. The relative 
timing between driver and probe pulse was varied using an optical delay line. 
The probe, after passing transversely through the plasma, was split and directed 
into shadowgraphy and interferometry channels. These gave information not 
only about plasma creation and evolution, but also provided the in-situ neutral 
density profile. Along the optical axis, ~ 0.7 mm above the nozzle edge, the 
density had an approximately triangular density profile, going from zero to 
maximum density over a length of ~ 825 /im. The electron (and thus proton) 
density could be varied up to a maximum of rie ~ 10^^ cm~^ 10 ncr- The 
ion beam was characterised with a magnetic spectrometer, which dispersed the 
protons by their deflection in a transverse magnetic fleld. The aperture of the 
spectrometer was a = 0.6 mm diameter pinhole. The dispersed protons were 
detected with a (polyvinyltoluene) scintillator screen which light emission was 
calibrated to the dose of energy deposited by protons impinging on it. The 
scintillator was reimaged on to a Andor EMCCD camera. 

0.3. Data unfolding 

Proton energy spectra were unfolded by first removing hot-spots (due to x- 
rays), background subtracting, and then integrating vertically to give a line-out. 
An absolute residuals optimisation was used to find the trial spectra which best 
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reproduced the measured signal after convolving with the instrument function, 
which was taken from the original vertical spread of the signal. We note that 
due to the small acceptance angle of the spectrometer (9.8 x 10~^ sr), the 
transport through the spectrometer has little effect on the signal spread. Density 
information was obtained from the interferograms, by first obtaining the phase 
with reference to a background image. An Abel transform was then used to 
obtain the density profile from the phase change, assuming cylindrical symmetry 
around the laser axis. 
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